A Relativistic Separable Potential to Describe Pairing in Nuclear 

Matter 



B. Funke Haas, T. Frederico, B. V. Carlson 

Departamento de Fisica, Instituto Tecnologico da Aerondutica - CTA, 
12.228-900 Sao Jose dos Campos, Sao Paulo, Brazil 

F.B. Guimaraes 
EAN - Instituto de Estudos Avangados - CTA, 
12.228-840 Sao Jose dos Campos, Sao Paulo, Brazil 



Abstract 

Using the Dirac-Hartree-Fock-Bogoliubov approximation to study nuclear 
pairing, we have found the short-range correlations of the Dirac 1 So pairing 
fields to be essentially identical to those of the two-nucleon virtual state at all 
values of the baryon density. We make use of this fact to develop a relativistic 
separable potential that correctly describes the pairing fields. 
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I. INTRODUCTION 



We have recently studied nuclear pairing in meson-exchange models of the nuclear in- 
teraction using a Dirac-Hartree-Fock-Bogoliubov (DHFB) approximation to nuclear matter 
IJTIIJ]. An important conclusion of this study is that the short-range 1 So pairing correlations 
in nuclear matter are essentially identical to the short-range correlations of the two-nucleon 
virtual state, due to the dominance of the virtual state in the 1 So channel. In non-relativistic 
calculations, Khodel, Khodel and Clark || also noted a close similarity between the pairing 
gap function and the short-range vertex function of the virtual state M, which prompted 
them to suggest that the virtual state be used as the starting point in calculating the gap 
function. In our relativistic calculations, however, the similarity of the Dirac pairing fields 
and the vertex functions is so striking that we can claim that the two are essentially identical. 
We make use of this quasi- identity to develop a separable approximation to the relativistic 
pairing potential that correctly takes into account the high momentum contributions of the 
short-range two-nucleon correlations. 

Another indication that a separable potential can describe nuclear pairing correlations 
comes from the usual assumption of a constant pairing gap in finite nucleus calculations. The 
analysis of the HFB equations in the zero-range limit of the two-nucleon interaction shows 
that the separable potential approximation becomes an exact formulation of the pairing 
problem in the zero-range limit as the pairing potential then becomes independent of the 
baryon momentum. The pairing field also become independent of the baryon momentum 
and vary only with respect to the baryon density of the nuclear matter, which is equivalent 
to a constant pairing gap assumption. Thus, the separable potential approximation for the 
pairing field can be seen as the simplest extension, to finite-range two-nucleon interactions 
and all baryon densities, of the zero-range limit of the HFB-equations. It gives the simplest 
improvement to the usual assumption of a constant pairing gap. 

The aim of this paper is to present a rank-one, separable interaction that accurately de- 
scribes the DHFB pairing field at all densities. In section II, we clarify the close relationship 
between the pairing field and the two-nucleon 1 So virtual state vertex function and use this 
relationship to define the separable interaction. In Section III, we demonstrate the goodness 
of the separable interaction through comparisons with exact DHFB calculations of nuclear 
matter. 
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II. THE FORMALISM 



The bound-state correlations in a two-particle system can be roughly classified as asymp- 
totic ones and short-range ones. The asymptotic correlations are determined principally by 
the binding energy, while the short-range ones depend on the high-momentum components 
of the wave function. This can be seen by examining the manner in which a bound pair 
appears in the two-body T-matrix, T(E). The T- matrix satisfies the integral equation 

T(E) = V + VG (E)T(E), (1) 

where V is the two-body interaction and Gq(E) the free two-body propagator. A bound 
state corresponds to a pole in the T-matrix at a negative value of the energy, — e&. We can 
decompose the T-matrix, in this CclSC, clS 

T(E)=T Lrt+T c (£), (2) 

— e b — t, 

where T is the bound-state vertex function and T C (E) is the continuum component of the 
T-matrix (and, possibly, the contributions of other bound state poles). Substituting the 
decomposition in the integral equation, Eq. (U), we verify that the vertex function satisfies 
the equation 

T = VG (-e b )T. (3) 

The vertex function of the two-nucleon bound state is closely related to its wave function, 
if} = Go(—€b)T, which satisfies the differential equation 

(-66 - H )tP = (-e 6 - H )G (-e b )T — V — VG (-e b )T = Vtfj , (4) 

where Hq is the Hamiltonian of the free two-nucleon system. 

Here we see the rough division of the two-nucleon correlations in the vacuum. The 
asymptotic correlations are determined by the singularity of the Green's function Gq(— e b ) 
while the short-range ones are contained in the high momentum components of the vertex 
function, T. This analysis can also be extended to the nucleon pairing correlations in nuclear 
matter. One can directly associate the bound-state vertex function T of the T-matrix with 
the pairing field of the HFB approximation, A, and the adjoint vertex function, V* , with 
the time-reverse conjugate pairing field, A. To demonstrate this, we develop the relation 
between T and A in a formal way, by comparing the equations that define the two. (In a 
similar manner, the pair wave function ip can be associated with the residue of the anomalous 
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propagator F in the complex energy plane 0|| . However, we will not develop this association 
here.) 

In the Dirac-Hartree-Fock-Bogoliubov (DHFB) approximation to pairing in nuclear mat- 
ter, the self-consistency equation for the pairing field can be written as [|],|| 

A ( k ) = i E / -0ji A ^Df(k - q)F(q)BAj^ , (5) 

where, Aj a represents the meson-baryon vertex of meson j (for the a meson, A aa = ig s 1, 
for the uj meson, A^ a = —igv^a-, etc.), Dj 13 represents the propagator of meson j and 
the Greek letters represent any necessary Lorentz and/or isospin indices. The matrix B 
relates transposed quantities to the complex conjugates of time-reversed ones. It is given 
by B = T2 <8> 75C, where the Pauli matrix T2 acts in the isospin space and C is the charge 
conjugation matrix. 

The anomalous propagator F(q) in the self-consistency equation, Eq. (|3p, is one of the 
number-non-conserving components of the baryon propagator Sp(q), which, in the Gorkov 
formalism, takes the form 

WH c,,)fW |, (6) 

' F(q) G(q) 1 

in which G(q) is the usual number-conserving propagator and F(q) and G(q) are the corre- 
sponding time- reversed propagators. The baryon propagator Sp(q) can be written in terms 
of the self-energy and pairing fields, S(g) and A (q), respectively, as 

Sr( , )= i^- M : E ('> + « A <'> V 1 , (7 ) 

A(,) A + M + Y.(q)-n) 

where \i is a Lagrange multiplier used to fix the average baryon density. 

The form of the 1 So pairing field can be determined from the hermiticity and antisymme- 
try properties of the Lagrangian density, as well as the translational, rotational and isospin 
symmetries of symmetric nuclear matter, as is shown in Ref . [|IJ . The general form obtained 
for the field is 

A(k) = (A s {k) - 70 A (fc) - 7o7 ■ k A T {k))r -h. (8) 

where k = (ko,k) is the baryon 4-momentum, k = k/\k\ and the isospin orientation n is 
arbitrary. We note that f ■ h = T2 corresponds to the standard case of proton-proton and 
neutron-neutron pairing. 
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We can decompose the self-consistency equation for the pairing field, Eq. ([|), into equa- 
tions for the components of Eq. (|8]). After integrating over energy and angle and neglecting 
retardation effects, the component equations take the form, 

A s (k, k F ) = -L / A q 2 dqV s (k, q)F s (q, kp) , 

A (fc, kp) = / q 2 dqV (k, q)F (q, k F ) , (9) 

Z7T Z Jo 

A T (k, k F ) = [ A q 2 dqV T (k, q)F T (q, k F ) , 
2ir z Jo 

where k = \k\ and g = \q\. The pairing potentials, Vg, V and (given in the Appendix), 
are functions of the coupling constants and meson masses, while A is a cutoff in the baryon 
momentum. The Fermi momentum kp is defined through its conventional relation to the 
baryon density, that is, pB = jkp/Sii 2 , where 7 is 2 for nuclear matter and 1 for neutron 
matter. The components of the anomalous propagator are defined as 

F 3 (k) = ±Tr[f-hF(k)), 

F {k) = ^Tr[ l0 T-hF{k)), (10) 

F T (k) = -\Ti[ 1 ^-kr-nF(k)}. 
o 

These are evaluated in the Appendix in the case in which their negative-energy Dirac sea 
contributions are neglected. 

The components of the anomalous propagator, Fg, F Q and Fp carry the information of 
the density-dependent nucleon self-energy and pairing mean-fields. The pairing potentials 
Vs, Vq and Vp also possess a density dependence associated with the retardation terms, which 
account for the finite meson propagation velocity. In Eqs.(|^), these retardation terms are 
neglected since their effects on the pairing fields have been found to be small Q. The den- 
sity dependence remaining in the self-consistent pairing equations, Eqs.@, is thus entirely 
contained in the components of the anomalous propagator. 

When the baryon density ps tends to zero, both the self-energy T,(q) and the pairing 
field A(g) also tend to zero. The anomalous propagator F(q) can then be well approximated 
by 

F ^)^- i M ] — -~ A (g) t- M 1 -, PB^O, (11) 

fa — M + [ijq + it] jq + M - fir/ Q - vq 

an approximation that becomes exact at zero density. The vacuum limit of the pairing field 
thus satisfies the homogeneous equation 
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A ( k ) = *E / 7^4 A ^f (* - g) ■/ fl/f .- A (g) ^IT--^ 51 • ^ 

^ 7 (27r) 4 j /i7o+ fy - M + ir] A*7q— jq - M + ir) Ji 

(12) 

Turning now to the Bethe-Salpeter equation, we write the bound-state vertex function 
as r(fc, P), where P is the center-of-mass four momentum of the baryon pair and k is their 
relative four momentum. The Bethe-Salpeter equation for the vertex function, in the ladder 
approximation, is often written as 

T(k,P) J ^A lja AwDf(k-q)G 01 (P/2 + q)G 02 (P/2-q)T(q,P), (13) 

where the free baryon propagator Go(q) is given by 

w-T^m- (14) 



We can rewrite this in a matrix form, making use of the matrix B, as 
T(k, P)P f = < E / 4^i A ^ D f( k " <?) G o( P / 2 + 9) r fo> P)B^G (P/2 - q)BAj^ . (15) 



Comparing this expression with the self-consistency equation for the vacuum pairing field, 
Eq. fljjD, we find the two are identical if we 

• evaluate the vertex function in the center-of-mass frame, taking 

P Q = 2fi5 a0 , (16) 

where ix is the Lagrange multiplier of the Gorkov propagator Spi^q) of Eq. (|7|), and 

• associate the pairing field and the vertex function as 

A{k) = T{k)B\ (17) 

where we now suppress the center-of-mass dependence of the vertex function. 

We thus conclude that the self-consistency equation for the vacuum pairing field is identical 
to the ladder approximation to the Bethe-Salpeter equation in the center-of-mass frame. 

An analysis of the symmetry and antisymmetry properties that we expect of the center- 
of-mass Bethe-Salpeter vertex function, T(k), leads us to the following form, 

r(fc)pt = f Ts (k) - l0 r (k) - 70 7 • kv T {k)) f-n, (is) 
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which is identical to that of the pairing field A(fc) given in Eq. fl8|). We substitute this 
vertex function in the Bethe-Salpeter equation, Eq. (|T5|), integrate over energy and angle, 
and neglect retardation effects and the contribution of the negative energy states, just as we 
have done in the self-consistency equation of the pairing field. We find, 

M ; Jo M ,q Un 2 E a (M B ~2E q ) y ' 



T (k) = I Vo (k, q )— 2 M B E q (M B - 2E q ) dq ' (19) 

v (V\ Pvtir ,Q 2 k q(M B T s -2MT ) 
T T (k) = -J o Vr(k, q )^ MBEq{MB _ 2Eq) d q , 

where k = \k\ , q = \q\, A is a cutoff in the baryon momentum, and E q = y/q 2 + M 2 . The 
potentials V$, V , and V? are, of course, the same as those of the self-consistency equations 
for the pairing field, Eqs. and are given in the Appendix. The center-of-mass energy of 
the bound two-nucleon pair, that is, the bound two-nucleon mass, M B , is given by 

M B = 2fi = 2M - e B . 

In the Appendix, the equivalence between the vacuum pairing equation and the ladder 
approximation to the Bethe-Salpeter equation is once again demonstrated, by reducing the 
explicit form of the components of the pairing equation, Eqs. to the components of the 
Bethe-Salpeter equation, Eqs. flT9p, in the vacuum limit. 

The above expressions are valid for a bound two-body system in the vacuum. If the pole of 
the T- matrix corresponds to a virtual state (anti-bound state), as in the case of the physical 
1 So two-nucleon channel, one must redefine the integrands through an analytic continuation 
into the second sheet of the complex energy plane . The analytical continuation introduces 
an additional term, which takes into account the discontinuity of the propagator across the 
cut in the energy plane. In the components of the Bethe-Salpeter equation, Eqs. (|19D , in 
which the negative-energy states have been neglected, the analytically continued equations 
can be obtained through the substitution, 

' ' - 2ixi5(M v - 2E q ) , (20) 



M B - 2E q M v - 2E q 
where M v is the mass of the two-nucleon virtual state 



M v = 2M — e v = ^M 2 -\q v \ 2 

with q v = —i\q v \. 
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The equations for the components of the vertex-function then become 
T s (k) = T s (k) v - {V s (k, q) (E q T s + MT )} q=qv , 

T (k) = T (k) v + — \ V (k, q) — ) (21) 



TT 



q=qv 



r T (fc) = T T (k) v + — <j kV T (k, q) — ) (22) 
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/ <j=<j„ 

The functions T s (k) v , T (k) v and Y T {k) v are given by Eqs. (|T9"D with the mass of the bound 
state Mb substituted by that of the virtual state M v . 

As stated in the introduction, we have found that the momentum dependence of the 
components of the pairing field A, calculated with various sets of interaction parameters, 
shows a very small variation over a wide range of values of the baryon density and a wide 
range of cut-offs (A between 2.5 fm _1 and 15.0 fm _1 ) 0. As an example, in Fig. [TJ we 
show the momentum dependence of the principal components of the pairing field, Ag(k) 
and A (k), for several values of the Fermi momentum, which were obtained using the a — uj 
parameters of Ref. with a momentum cutoff of A = 10 fm^ 1 . Note that the curves 
denoted by Fermi momentum zero correspond to the components of the vertex function of 
the x So virtual state in the vacuum. The small variation in the components as a function 
of the Fermi momentum suggests the use of the momentum dependence at one value of 
the density. Use of the vacuum momentum dependence is particularly convenient, since 
in this case a direct correspondence exists between the self-consistent pairing field and the 
Bethe-Salpeter vertex function. A direct comparison between the numerical solutions for 
the vertex function, T(k), and for the pairing field, A(k), gives us the following very good 
approximate relations between the components of the two fields, 

A s (k, k F )^d(k F ) T s (k), 

A (k,k F )^d(k F )V (k), (23) 
A T (k,k F )^d(k F ) F T (k), 

where d(k F ) is a density- dependent factor determining the overall magnitude of the pairing 
fields. 

We emphasize that, in the nonrelativistic approach of Khodel, Khodel and Clark, the 
virtual state vertex function serves only as a good starting point for calculating the pairing 
field , while in our relativistic approach, the two are essentially identical in form. We can 
explain the difference between the two approaches by analyzing the effective gap function 
for the positive energy DHFB states, which takes the form 
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where E* = V k* 2 + M* 2 and the effective momentum and mass of the quasi-nucleons, k* 
and M*, are defined in the Appendix. The gap function A 9 (fc) is the quantity whose role is 
closest to that of the non-relativistic pairing field Jl|. Like the non-relativistic self-energy, 
the gap function is the difference between two larger relativistic quantities, A and A s . 
(Ay is two to three orders of magnitude smaller than the other components.) The effective 
mass M* that appears in the gap function varies strongly with the nuclear matter density. 
Although the form of the components of the Dirac pairing field remain almost constant with 
density, the gap function changes, due to its dependence on M* . The density independence 
of the form of the relativistic pairing field thus does not carry over to the nonrelativistic 
field. 

The invariant momentum dependence of the components of the pairing fields suggests a 
separable form for an effective potential, in which each of the components of the interaction, 
Vs, V and Vt, is proportional to the corresponding component of the vertex function, T s , 
To, and IV . Such potentials are well known from the nonrelativistic treatment of the two- 
nucleon problem. It is such a separable form of the non-relativistic pairing potential that 
Khodel, Khodel and Clark suggest as a starting point for calculating the nonrelativistic gap 
function ||. Here, we use a form appropriate for the components of the relativistic pairing 
field that also preserves the symmetry in the two momenta, k and q, of the potential, taking 

Vs(k, q ) = ^T s (k)T s ( q ) , 

V (k iq ) = ^T (k)T (q), (25) 
kV T (k,q) = ^-T T (k)T T (q) , 

where the potential strengths, A5, Ao, and At, are obtained through the simultaneous solution 
of Eqs. (|1^) and (|25"D , in the case of a bound state, or Eqs. fl2~T|) and (^) , in the case of a virtual 
state. We note that, although we have included the equation for the tensor component IY 
in Eqs. ( p3| ) and (^5|), we have neglected this component in our numerical calculations. The 
inclusion of the tensor term offers little further calculational difficulty but, being extremely 
small, has little effect on the rest of the calculation. 

For any given meson-exchange potential, we have shown in Ref. @] that the size of the 
pairing gap, for low baryon densities, is correlated with the energy of the virtual state of 
the T- matrix, e v . The value of e v , in turn is related to the baryon momentum cut-off, A, in 
the integrals of the self-consistency equations. Consequently, fixing e v at its physical value 
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in the Bethe-Salpeter equation for T determines a (usually unique) value for A. 

The separable potential introduces two additional parameters, Xs and A , when the tensor 
term is neglected. These parameters, together with the momentum cutoff, A, are fixed by 
requiring that the energy of the two nucleon state and the ratio of the vertex functions be 
equal to the respective ones obtained by solving the Bethe-Salpeter equation with a meson- 
exchange potential. Together, the parameters A^, Ao and A introduce an ambiguity in the 
determination of the separable potential. Our calculations have shown that, within a wide 
range of values of the cutoff A, we can find values of As (A) and Ao(A) that fix the virtual 
state energy and the vertex functions at their physical values. We regard this as a positive 
feature of the separable potential, as it permits us to fix one of the parameters, usually A, 
at an arbitrary, and possibly convenient, value, and then solve for the others. 

The only task remaining to complete the definition of the separable pairing potential is 
to extend the potential, which we have defined in the vacuum, to finite values of the baryon 
density. This is done by observing that, as the retardation terms in the pairing potentials 
have been neglected, the dependence of the pairing field on the baryon density is completely 
determined by the density dependence of the components of the anomalous density, Fg, F 
and Ft, in Eq.(||). Consequently, the separable pairing potential for finite baryon densities 
must be the same as the potential in the vacuum. The self-consistency equations for the 
components of the pairing field are thus given by Eqs. (|^) and (|25|), with the parameters 
Xs, A , and A (and, eventually, A^ as well) fixed by the vacuum solution. 

III. RESULTS AND CONCLUSION 

We have compared calculations of the pairing fields obtained using a separable interac- 
tion with the usual DHFB results for various sets of a — uj interaction parameters. The 
relativistic separable potential describes well both the gap parameter, A g , and the two 
large components, A and A s , at all baryon densities. It also reproduces the behavior of 
the pairing field A for increasing values of the momentum cut-off A in the self-consistency 
equations. 

In Fig. ^|we show the particular results for the cr — uj interaction of Bouyssy et.al. 0, for 
a cut-off of A = 3.8 fm -1 , obtained by adjusting the virtual state energy to its physical value. 
This potential has been chosen because its vacuum gap function, which is the nonrelativistic 
two-nucleon vertex function, is in good agreement with the nonrelativistic two-nucleon vertex 
function calculated with the Bonn-B potential ||. We note a difference of at most a few 
percent between the components of the relativistic fields obtained in exact and separable 
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calculations. The differences in the two gap functions are even smaller. Similar results are 
obtained for other sets of a — u interaction parameters [0 . 

As discusses above, it is possible to adjust the strengths As and A of the separable 
interaction to fit the experimental value of the two-nucleon virtual state energy as a function 
of the cutoff A. For each value of the A, the resulting values of the strengths As and Ao 
yield a separable potential which carries the physical information necessary for the correct 
evaluation of the gap parameter at all values of the density. Since the properties of the 
virtual state dominate the physics of 1 S pairing, we would expect the pairing fields that 
result to be independent of the cutoff used. We can verify the extent to which this is true by 
comparing the results obtained for d(kp) of Eq. (|23|), the common factor determining the 
overall magnitude of the pairing fields, at different values of the momentum cutoff A. We do 
this in Fig. || where we display the factor d{kp) as a function of the Fermi momentum for 
several values of the cutoff A. We find the form of the factor d(kf) to be fairly constant but 
its magnitude to decrease between 5 and 10 % as the momentum cutoff A increases from 3.3 
fm _1 to 10 fm _1 . 

In Ref . [J2J , a strong correlation between the magnitude of the pairing gap and the energy 
of the two-nucleon virtual state was observed in all the meson potentials studied there. The 
correlation is reproduced, for all values of the cut-off, when using the corresponding separable 
potential. In Fig. [|, we compare the correlation for the gap parameter, A g , obtained with 
the separable interaction with the exact result, using the o — uo parameters of Ref. [[j]]. The 
agreement between the two is excellent, confirming once again the goodness of the separable 
approximation to the pairing interaction. 

We thus conclude that a separable interaction obtained by exploiting the close rela- 
tionship between the Dirac pairing field and the virtual state vertex function furnishes an 
excellent description of the Dirac pairing fields at all values of the baryon density. We next 
plan to implement a configuration-space version of the separable interaction in a DHFB code 
for finite nuclei calculations. We expect that such a potential will describe the short-range 
pairing correlations with greater accuracy than those obtained with the usual constant gap 
assumption. 
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APPENDIX: THE ANOMALOUS PROPAGATOR AND THE PAIRING 

POTENTIAL 

The functions F$, F and Ft in Eq.(6) are given by [|l|], 

F s {k) = {A s [-a + (3 + Ef + M* 2 + k* 2 ] + AJ(A§ - A 2 + A 2 T k 2 ) 

-2M*£* A - 2k*E%kA T } /(4/3w_) , 

F (k) = {A (a -(3-Ef- M* 2 + k* 2 ) + A*(A| - A^ + A 2 T k 2 ) 

+2M*E*A s -2k*M*kA T }/(4Puj , (Al) 



and 



F T (k) = {-kA T (-a + (3 + Ef - M* 2 + k* 2 ) - kA* T (A 2 s -A 2 + A 2 T k 2 ^ 

+2k*E^A s - 2M*k*A } /(4/9u;_) , 

where the asterisk indicates complex conjugation and k = \k\. The quantity c<j_ is 

u _ = y/a-p (A2) 

where 
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a = Ef + k* 2 + M* 2 + \A S \ 2 + |A | 2 + |A T | 2 

A = -Im(A T A*) , 
fa = EpM* + Re(A 5 A*) , 
(3 3 = Epk* + Re(A s A* T ) , 
04, = \M*A T + k*A \ . 



(A3) 



The effective momentum and mass of the baryons, k* and M*, are given in terms of the 
components of the nucleon self-energy, 



E(k) = E s (k)- l0 j: (k)-i-kj: v (k), 



(A4) 



as, 



k* = (1 + E v (k))k and M* = M + E s (fc) 



(A5) 



The Fermi energy is 



E F = Z (k)+n, 



where \i is the chemical potential. 

The pairing potentials that take into account the exchange of the mesons a and uo are 
given by 



V T (k,q) = 

where the functions 6 a and <p a are 



2k 2 q 



~^9% + \A 

\g% + \g 2 J u 



gl<t>* 



= IH a + 2k q \ 
\H a -2kq 



with 



(A6) 



(AT) 



H a = (u4k) - u.(q)Y - k 2 - q 2 - m 2 a , 



(A8) 



and 



Akq a ' 



(A9) 
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Here a indicates the mesons of the model, with a = a, 00. 
In the limit of zero baryon density, we have 



k* -> k, M* — > M, and ji — > E* F — > M s /2 , 

so that 

a->^ + ££, (3^M B E kl and w_ ^ ~ |M B - 2£ fe | , 

where = \//c 2 + M 2 . The components of the pairing field tend to zero. The zero baryon 
density limit of the pairing self-consistency equations are then 

q 2 E q A s -MA c 



a fiA I V/U \ V A 
As(k)=J Q Vs(k, q )— 2 Eq{MB _ 2Eq) d q 



A (k)=J o V (k, q )— 2 M B E q (M B - 2E q ) dq > 

V; Jo An 2 M B E q [M B — 2E q ) y ' 

which have exactly the same form as the equations for the components of the vertex function 
T given in Eqs.(|l9D. 
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FIGURES 




k (fm _1 ) 

FIG. 1. The momentum dependence of the two large components of the Dirac pairing field, A^ 
and Ao, are shown at several values of the nuclear matter density, using the parameters of Ref. 
[7]. The curve labeled as zero density corresponds to the vertex function of the virtual state. 
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FIG. 2. The pairing fields obtained in exact and separable DHFB calculations, using the a — u> 
parameters of Ref. [7], are shown. The two large components of the Dirac pairing field, Ag and 
Aq, and the gap parameter, A g , are shown as functions of the baryon density, for A = 3.8 fm -1 . 
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FIG. 3. The factor d(/cp) determining the overall magnitude of the pairing fields, calculated 
using the a — u> parameters of Ref. [7], is shown for several values of the momentum cutoff A. 
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FIG. 4. The magnitude of the pairing gap as a function of the energy-equivalent momentum 
of the two-nucleon virtual state from exact and separable DHFB calculations, using the a — u 
parameters of Ref. [7], are shown. The positive and negative values of the abscissa correspond to 
bound and virtual two-nucleon states, respectively. 
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